A Large Atom Number Metastable Helium Bose-Einstein Condensate 
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We have produced a Bose-Einstein condensate of metastable helium ( 4 He*) containing over 1.5 x 10 7 
atoms, which is a factor of 25 higher than previously achieved. The improved starting conditions 
for evaporative cooling are obtained by applying one-dimensional Doppler cooling inside a magnetic 
trap. The same technique is successfully used to cool the spin-polarized fermionic isotope ( 3 He*), 
for which thermalizing collisions are highly suppressed. Our detection techniques include absorption 
imaging, time-of-flight measurements on a microchannel plate detector and ion counting to monitor 
the formation and decay of the condensate. 



Ten years after the first experimental realization of Bose- 
Einstein condensation (BEC) in dilute, weakly interact- 
ing atomic systems the field of degenerate quantum 
gases has developed into a major area of research. For 
most elements it has not yet been possible to produce 
Bose-Einstein condensates containing large numbers of 
atoms. Only for hydrogen, sodium and rubidium have 
condensates with more than 10 6 atoms been realized 0. 
Large condensates provide a better signal-to-noise ra- 
tio, allow a study of both the collisionless and the hy- 
drodynamic regime and are especially useful for sym- 
pathetic cooling and atom laser applications. In this 
realm metastable atoms are of particular interest, of- 
fering alternative detection methods due to their high 
internal energy. Absorption imaging is the technique 
most frequently applied to detect and measure a BEC. 
Metastable helium (2 3 Si state) is the only species for 
which detection of the condensate has been performed 
using a microchannel plate (MCP) detector 0; the same 
detector was also used to measure ions produced by 
Penning ionization (while absorption imaging was not 
used). In a recent experiment, after pioneering work 
with metastable neon [4|], Schellekens et al. used a po- 
sition sensitive MCP detector to observe the Hanbury 
Brown and Twiss effect both above and below the BEC 
threshold In the second experiment in which BEC 
was realized with He* 0, optical detection was used 
and up to 6 x 10 5 atoms could be condensed. Recently, 
this group reported a high precision measurement of the 
metastable helium scattering length a, which was per- 
formed using a two-photon dark resonance • The value 
a = 7.512 ± 0.005 nm is favorable for experiments with 
the fermionic isotope 3 He*. It ensures a stable ultracold 
4 He*/ 3 He* boson- fermion mixture, as the inter- isotope 
scattering length will be large and positive ||. Large 
numbers of 4 He* atoms all along to the critical temper- 
ature provide an efficient reservoir for sympathetic cool- 
ing and will facilitate the production of degenerate 3 He* 
clouds with large numbers of atoms. 

In this letter we present an experiment that com- 
bines the various detection methods used previously 0,0 
and describe the realization of a BEC of 4 He* con- 
taining more than 1.5 x 10 7 atoms. This large im- 



provement is primarily due to the application of one- 
dimensional Doppler cooling inside the magnetic trap 
rather than three-dimensional Doppler cooling prior to 
magnetic trapping. Doppler cooling of polarized atoms 
was originally proposed for atomic hydrogen [t| , and has 
recently been demonstrated for optically dense samples of 
magnetically trapped chromium atoms 10]. Compared 
to the laser cooling methods we investigated previously 
1 1 1 L Il2j|. this configuration is more efficient and simple. 

The experimental setup is an extended and improved 
version of our previous setup In short, we start 

with a beam of metastable atoms produced by a liquid- 
nitrogen cooled dc-discharge source. The atomic beam is 
collimated, deflected and slowed by applying laser beams 
resonant with the 2 3 Si — ► 2 3 P2 transition at 1083 nm. 
Typically 2 x 10 9 atoms are loaded into a magneto-optical 
trap (MOT) at a temperature of 1 mK. Since our previ- 
ous experiments |l lj we have installed a new ultrahigh 
vacuum (UHV) chamber and magnetic trap. The coils for 
our cloverleaf magnetic trap are placed in water-cooled 
plastic boxes and positioned in re-entrant windows. In- 
side the UHV chamber two MCP detectors (Hamamatsu 
F4655) and an rf-coil are mounted. The first MCP de- 
tector is positioned ~10 cm from the trap center and 
attracts positively charged ions produced in Penning ion- 
izing collisions: He* + He* — > Hc+ + Hc(l *S) + e~ (or 
He* + He* — > He2 + + e _ ). These ionization processes are 
the primary loss mechanisms in cold clouds of He*. A 
second (identical) MCP detector shielded by a grounded 
grid is positioned 17 cm below the trap center and de- 
tects neutral He* atoms that fall upon it. This detector is 
mounted on a translation stage and can be displaced hor- 
izontally to allow a vertical laser beam to pass through 
the trap center for absorption imaging. Absorption imag- 
ing of the MOT cloud determines the number of atoms in 
the MOT with an accuracy of about 20% 0] ; this is used 
to calibrate the He* MCP detector. When the MOT is 
loaded we switch off all currents and laser beams. In an 
applied weak magnetic field we spin-polarize the cloud 
(atoms are pumped into the m — +1 magnetic sublevel) 
and switch on the currents of the cloverleaf magnetic 
trap. Typically ^60% of the atoms is transfered from the 
MOT into the magnetic trap in this procedure. We oper- 
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FIG. 1: Time-of-flight signals of 4 He* atoms released from 
the magnetic trap, with and without one-dimensional Doppler 
cooling. The apparent signal increase after Doppler cooling is 
due to the increased fraction of atoms that is detected at lower 
temperature. The line is a fit assuming a Maxwell-Boltzmann 
velocity distribution. 



ate the cloverleaf trap at a bias magnetic field Bq—24 G 
to suppress excitation of depolarizing transitions in the 
subsequent one-dimensional Doppler cooling stage. 

One-dimensional Doppler cooling starts at the same 
time as the cloverleaf trap is switched on. It is im- 
plemented by retroreflecting a weak circularly polarized 
laser beam along the (horizontal) symmetry axis of the 
magnetic field. During the cooling pulse the tempera- 
ture decreases reducing the size and increasing the optical 
thickness. Cooling in the radial directions relies on reab- 
sorption of spontaneously emitted red-detuned photons 
by the optically thick cloud. Other possible energy redis- 
tribution mechanisms are collisional thermalization and 
anharmonic mixing. While the collision rate increases 
from 1.5 s _1 to 20 s _1 during Doppler cooling, anhar- 
monic mixing is negligible in our trap. With a laser de- 
tuning of one natural linewidth below the resonance fre- 
quency for an atom at rest in the center of the trap and 
an intensity of 10 _3 / sa t (I S at = 0.17 mW/cm 2 ), optimum 
cooling is realized in 2 seconds. In a separate experiment 
10 8 3 He* atoms were loaded into the magnetic trap and 
cooled using the same technique. For identical fermions 
s-wave collisions are forbidden, while the contribution of 
the higher-order partial waves is highly suppressed in this 
temperature range for He*. We observed a temperature 
decrease from 1 mK to 0.15 mK, which suggests that re- 
absorption of red-detuned photons scattered by atoms in 
the cloud is the main cooling mechanism in the radial 
direction. 

Figure ^ shows two TOF traces, illustrating the effect 
of one-dimensional Doppler cooling in our cloverleaf mag- 
netic trap. We typically trap N = 1.5 X 10 9 atoms, which 
are cooled to a temperature T = 0.13 mK, three times 
the Doppler limit. This implies a temperature reduction 
by a factor of 8 and an increase in phase-space density by 
a factor of ^600, while practically no atoms are lost from 



the trap. For comparison, reaching this temperature by 
means of rf-induced evaporative cooling would result in 
the loss of ~90% of the atoms from the trap. In pre- 
vious experiments [ill Il2j| we applied three-dimensional 
Doppler cooling or a two-color magneto-optical trap to 
improve the starting conditions for evaporative cooling. 
One-dimensional Doppler cooling provides lower temper- 
atures, higher phase-space density and is easier to im- 
plement. At this point the lifetime of the atoms in the 
magnetic trap is about 3 minutes, limited by collisions 
with background gas. To compress the cloud we reduce 
the bias field in the trap center to 3 G in 2.5 seconds, 
which increases the temperature to 0.2 mK. The param- 
eters of our magnetic trap then are modest: the axial 
and radial trap frequencies are lu z /2tt = 47 ± 1 Hz and 
oj±/2it = 237 ±4 Hz respectively. Higher frequencies are 
possible but are not required to achieve BEC. Our pro- 
cedure is similar to previous experiments on evaporative 
cooling in our group fl3| . After compression we cool the 
gas by rf-induced evaporative cooling in 12 s to BEC, 
which is achieved at a temperature of ~ 2/iK. We apply 
a single exponential rf ramp, starting at 50 MHz. The 
frequency decreases to zero but the ramp is terminated 
at around 8.4 MHz. Shorter ramps, down to 2 s, also 
produce a condensate, albeit with fewer atoms. 

The most sensitive method to detect BEC is TOF anal- 
ysis of the expanding cloud on the He* MCP detector. A 
typical TOF signal, obtained in a single shot, is shown 
in Fig. The double structure is a combination of 
the inverted parabola indicating the presence of a BEC 
released from the harmonic trap and a broad thermal dis- 
tribution. This signal is used to determine the number of 
atoms in the condensate as well as in the thermal cloud. 
In contrast to the Orsay experiments 0- Q- all atoms 
stay in the m = +1 sublevel during the trap switch-off. 
Applying the MCP calibration, the area under the fit- 
ted curve determines the number of atoms that have hit 
the detector. When we consider a thermal cloud, this 
number is only a small fraction of the total number of 
thermal atoms Nth- Therefore the determination of Nth 
relies upon the measured temperature and the MCP cal- 
ibration. The condensate expansion, determined by the 
mean-field interaction energy, is much slower. Thus the 
condensate will fall completely on the detector sensitive 
area (diameter of 1.45 cm), allowing us to measure the 
number of condensed atoms No using the MCP calibra- 
tion alone. The maximum number of atoms in the con- 
densate deduced in this way is 1 x 10 7 . This number is 
an underestimate due to MCP saturation effects, which 
will be discussed below. By applying an additional mag- 
netic field pulse we compensate small residual field gra- 
dients, which otherwise lead to a slight deviation from 
free fall expansion. With stronger field pulses we can 
also push the cloud towards the detector and this way 
realize shorter expansion times. The model used to fit 
the time-of-flight signals of the partly condensed clouds 
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FIG. 2: Observation of BEC, (a) on the He* MCP detector; 
the dashed fit shows the condensed fraction and the dashed- 
dotted fit the broader thermal distribution, (b) on a CCD 
camera; after an expansion time of 19 ms a round thermal 
cloud surrounding a cigar-shaped condensate is visible, (c) 
on the ion MCP detector; the condensate starts to grow at 
t = —0.95 s, at t = —0.45 rf ramp ends and at t = the trap 
is switched off. 



(Fig. |2Jt) is a combination of a thermal part (Bose distri- 
bution) and a condensate part (parabolic distribution). 
The chemical potential fj,, the number of atoms and the 
temperature T are the free parameters of the fit; effects of 
interactions are not included in the function used for the 
thermal part. In the Thomas- Fermi limit the chem- 
ical potential is given by [ftl 2 — 15ti 2 m 1 / 2 2~ 5 / 2 N u> 3 a, 
where Ti is Planck's constant divided by 27r, u> is the ge- 
ometric mean of the trap frequencies, m is the mass of 
the helium atom and a — 7.512(5) nm is the scatter- 
ing length. A maximum value of /i extracted from the 
fit of the TOF signal is ~ 1.3 x 10~ 29 J, which corre- 
sponds to 5.1 x 10 7 atoms in the condensate. A possible 
cause for the discrepancy between the number of atoms 
determined from the integrated signal and from the mea- 



surement of the chemical potential may be saturation of 
the MCP when detecting a falling condensate (peak flux 
of ~ 10 9 atoms/second); this leads to an underestimation 
of No as well as fi. Another possible cause is distortion of 
the velocity distribution during the trap switch-off and 
influence of remaining magnetic field gradients on the ex- 
pansion of the cloud. This may lead to an overestimation 
of fi, and therefore also of Nq. 

When the MCP detector is shifted horizontally, we 
can detect the condensate on a CCD camera. A weak 
(I = 10 _1 / sat ), 50 [is long on-resonance laser pulse is 
applied to image the shadow of the atoms on the CCD 
camera for which a quantum efficiency of ~1.6% is mea- 
sured at 1083 nm. As expected, the condensate ex- 
pands anisotropically while the thermal cloud shows an 
isotropic expansion (Fig.JSJj). Absolute calibration of the 
number of atoms at ~ //K temperatures could not be per- 
formed by optical means. The analysis of the absorption 
images, taken between 1 ms and 70 ms after the trap 
was switched off, shows that the condensate expansion 
deviates from the theoretical predictions : it expands 
faster than expected in the radial direction and slower 
in the axial. From these observations we conclude that 
the expansion of the cloud is influenced by magnetic field 
gradients during the switch-off of the trap. A difference 
in the switch-off times of the axial and radial confinement 
could cause an additional imbalance in the redistribution 
of the condensate interaction energy between the two di- 
rections. This may influence the measurements of both 
the chemical potential and the temperature. In order to 
check if the interaction energy is conserved, we extract 
the asymptotic kinetic energy gained durin g th e expan- 
sion from absorption images of the cloud [lfj. In the 
Thomas-Fermi approximation this so-called release en- 
ergy should equal the interaction energy in the trap. We 
obtain iVo = 4 x 10 7 from this analysis assuming that no 
extra energy is added to (or taken from) the system dur- 
ing the trap switch-off. This is not exactly fulfilled in our 
case as switching is not fast enough to ensure diabaticity. 

To verify our TOF signal analysis, we plot the chem- 

2/5 

ical potential as a function of N using data obtained 
from the MCP measurements (here Nq is the number of 
condensed atoms measured by integrating the MCP cur- 
rent). The data points lie on a straight line which goes 
through zero with a slope larger than expected, meaning 
that either [i is overestimated or Nq is underestimated. 
The former is supported by the analysis of the absorp- 
tion images so we correct fi. The corrected data points 
as well as the theoretical line are presented in the inset 
of Fig. |3 The plot also shows that the MCP detector 
saturates, when the number of atoms in the condensate 
exceeds ~ 10 6 . When we now extract the number of 
atoms from the measured chemical potential, after the 
correction for the distortion during the trap switch-off, 
we find TYo = 1.5 x 10 7 in our largest condensates. This 
number is still a lower limit, as the analysis assumes that 
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FIG. 3: Decay of a quasi-pure BEC (circles) and a BEC in 
the presence of a large (Nth = No) thermal fraction (squares). 
The dashed curve represents the atomic transfer model |19|. 
with two- and three-body loss rate constants obtained from 
a fit (full curve) to the decay of the quasi-pure BEC. Data 
points that lie above No = 10 6 are corrected for the saturation 
effects. Inset: Chemical potential (fi) as a function of N^ 5 . 
The same data as for the quasi-pure BEC decay are used in 
the plot. The value of [i, is multiplied by a constant (0.61), to 
bring the data points on the theoretical line (see text). 

/i is not affected by saturation of the MCP detector. We, 
however, measure a reduction in \i when we push the 
BEC towards the detector, thus increasing the satura- 
tion problem. 

With a second MCP detector we observed the growth 
and decay of our condensate by counting the ions pro- 
duced during evaporative cooling. Due to the increase 
in density the ion signal increases, although the number 
of trapped atoms decreases. When BEC sets in, a sharp 
increase is expected, indicating the formation of a dense 
cloud in the centre of the trap Q . This is demonstrated 
in Fig. 0;, which shows the growth of the condensate as 
well as its decay. 

The dynamics of formation and decay of the conden- 
sate is an interesting aspect that was discussed and in- 
vestigated earlier to some extent |l8j . In our group a 
model was developed describing the deca y of the conden- 
sate in the presence of a thermal fraction . The model 
assumes thermalization to be fast compared to the rate 
of change in thermodynamic variables, so the system re- 
mains in thermal equilibrium during the decay. It was 
shown that under this assumption a transfer of atoms 
should occur from the condensate to the thermal cloud, 
enhancing the condensate decay rate. To verify this, we 
performed measurements of the BEC lifetime using the 
TOF signal. Due to the high detection efficiency it was 
possible to detect a condensate up to 75 s after it was pro- 
duced. Results of these measurements are summarised 
in Fig. |31 We fit the model to the experimental data for 
a quasi-pure BEC decay; two- and three-body loss rate 
constants are used as free parameters. Good agreement 
with the experiment is found for two- and three-body loss 



rate constants (3 — 2 x 10 -14 cm 3 s _1 and L — 9 x 10~ 27 
cm 6 s _1 respectively, which compares well with theory 
[l7j | . When we use the extracted values of (3 and L in our 
model for the decay of the condensate in the presence of 
a thermal fraction, the dashed curve included in Fig.[3]is 
obtained. The agreement with the experiment is good, 
so we can conclude that the model reproduces the data. 

To summarize, we have realized a condensate of 4 He* 
containing more than 1.5 x 10 7 atoms and studied its 
growth and decay by measuring the ion production rate 
in situ, observing its ballistic expansion by absorption 
imaging and by recording the time-of-flight signal on an 
MCP detector. The main ingredient that made this large 
atom number possible is one-dimensional Doppler cool- 
ing in the magnetic trap. We demonstrated that this 
technique can also be applied to cool spin-polarized he- 
lium fermions, where the Pauli principle forbids s-wave 
collisions. Combining both isotopes in one setup may al- 
low the observation of Fermi degeneracy in boson-fermion 
mixtures of metastable atoms. 
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